Majorana neutrino textures from numerical considerations: 
the CP conserving case 
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Phenomenological bounds on the neutrino mixing matrix U are used to determine numerically 
the allowed range of real elements (CP conserving case) for the symmetric neutrino mass matrix M^ 
(Majorana case). For this purpose an adaptive Monte Carlo generator has been used. Histograms 
are constructed to show which forms of the neutrino mass matrix M^ are possible and preferred. 
We confirm results found in the literature which are based on analytical calculations, though a few 
differences appear. These cases correspond to some textures with two zeros. The results show that 
actually both normal and inverted mass hierarchies are still possible at 3a confidence level. 
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I. INTRODUCTION 

The standard neutrino theory involves diagonalization of the neutrino mass matrix M^, by use of the mixing matrix 
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m^ag = U^M^U. 



(1) 



The matrices nidiag and U are restricted by experimental data. Merging information from neutrino oscillation 
experiments with neutrinolcss double beta and tritium decay data |l|-|j|, the absolute neutrino masses can be found 
to be in the range: 



TOi<2.2eV, |TOi-mj| < 0.05eV, i,j = l,2,3. 
while U is restricted by the global neutrino oscillation analysis Q , see Table Q] 
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TABLE I: The allowed absolute values of the neutrino mass squared differences Ain^a, Arriai and the allowed absolute values 



of the neutrino mixing matrix elements \Uij\ 



and Oi are the central values and the 3o" uncertainties, respectively. 



However, what is a possible structure of the basic neutrino mass matrix Af^? Its knowledge could definitely shed a 
light on the mechanism of neutrino mass generation and is a subject of intensive studies Q-llfl- We concentrate on 
theories where quark and lepton mixings are disconnected. Then a base with diagonal charged lepton mass matrix 
can be used and information on neutrino mixings is connected solely with neutrino mass matrix M^. There are many 
papers devoted to the problem of determination of M^^. Basically, they can be divided into two groups where different 
methods to solve the problem are used. The first, "top-down" method relies on theoretical considerations of possible 




texture zeros and global symmetries which seems to arise from the neutrino mass matrix structure, e.g. [T^, [lll|.[T3|. 
There are also papers where the relation Eq.^is used to determine M^ analytically in terms of neutrino masses and 
values of neutrino mixing matrix elements, e.g. [131 , [ill , 13 ' ' [i3 • The second, "bottom- up" method relies on 
numerical analysis. One can perform numerical diagonalization of many neutrino mass matrix textures leaving only 
these which are in agreement with present experimental data, e.g. jl^, [T^l, [Ifl]. In this work we follow the second 
way: numerical solutions of neutrino mass textures are performed which are based on Monte Carlo simulations |27| . 

We have scattered a huge number of possible texture zeros, our analysis confirms most of the results based on 
analytic approach: we did not find numerical solutions for textures with number of zeros n > 3, we have found seven 
two zero textures which are in the range of current experimental results. Other eight two zero textures gave no results. 
For some discrepancies, see the next section. 

In this paper the real symmetric 3x3 neutrino mass matrix is analyzed. It means that we assume directly the 
Majorana nature of neutrinos and that the investigation is restricted to the CP conserving case. Certainly, it may 
seem to be a simple case, however, we believe that our visual presentation of results and conclusions show that even 
this case is nontrivial. The case with CP phases will be analyzed in a separate work. 

A general neutrino mass matrix Mi, which we analyze has the following form: 

(3) 

It has six parameters. Our approach is the following. First we consider the general case where all six parameters 
in Eq. |3| can be nonzero. Taking into account data in Table Q] we found the allowed range of these parameters. 
Histograms (frequency spectrums) for neutrino mass matrix elements (both for normal and inverted mass hierarchies) 
have been also constructed. These histograms give a nice visualization of possible solutions for M^, and give a hint 
which regions of parameters are in favor (e.g. for texture zeros). Next, a detailed analysis of the M^, textures are 
made where some of the parameters are exactly zero. Finally, conclusions are given and the Appendix contents details 
on the Monte Carlo algorithm used in the paper. 

II. THE GENERAL CASE 

Let us perform first numerical fitting for the general case of the M^, neutrino mass matrix (all parameters non-zero). 
To find numerical solutions, the matrix Eq. Owith some randomly generated numerical elements is diagonalized and 
the eigensystem is found: eigenvalues (neutrino masses) and normalized eigenvectors (columns of the neutrino mixing 
matrix U). Using eigenvalues, neutrino mass squared differences can be constructed: Am^i and Am^2- Both Arn^i, 
ATO32 and elements of U (let us call them Xi points) are compared with experimental values written in Table Q] (a; J*^"*). 
The rejection rule for a considered texture is defined in a form: 

i^±:i^ < A- (4) 

If inequality Eq. ^ is fulfilled, then the parameters Xi are possible solutions of the problem. To find the best values 
of parameters, the whole procedure is more refined, we send a reader to the Appendix A where the whole algorithm 
based on Adaptive Monte Carlo method is described. 

The parameter a gives an information about decreasing of the assumed error ai. If a = 1, then experimental data 
from the Table [H are taken and results are obtained at 3a confidence level. If a > 1 then we assume that errors 
shrink, however, we still assume the central values a;^^"* of Table Q] In this way a informs us how much Xi solutions 
are focused around their central values (see section III) . 

Let us note that modulus of the neutrino mixing matrix elements \Uij\ are analyzed and it is not possible to 
distinguish rotation symmetries of angles 0^, as a consequence, the figures included here have symmetric regions of 
allowed parameters. 

The algorithm is able to distinguish the normal mass hierarchy from the inverted one by finding only such combina- 
tions of neutrino mass matrix parameters (textures) which imply neutrino masses mi < m2 ^ ms (normal hierarchy) 
and ?7i3 -C mi < m2 (inverted hierarchy). In this way thousands of sets of allowed parameters can be find for a given 
hierarchy, these are plotted in Fig. (21 To allow for easier visualization of results, three dimensional projections of 
regions on the basic planes are shown (gray colours). In Fig. histograms for all M^ mass matrix parameters are 
shown. One can notice that histograms for normal mass hierarchy with a, 6, c parameters have their maxima around 
zero. The opposite situation takes place for inverse mass hierarchy Fig.^ Here the zero elements are connected with 
d, e and / parameters. These solutions suggest that textures with some parameters set to zero should be analyzed 
with a greater care. It will be done in the next section. 



III. TEXTURE ZEROS 

Here we will assume some parameters of the matrix M^ to be zero, exactly. It appears, that taking into account 
present experimental data Tabled we got in most cases solutions with both normal and inverse mass hierarchies. For 
an interesting case when a = (it means that the effective neutrinoless double beta decay mass parameter < m^^oi^ > ■• 
see e.g. [lj-|3 is zero for this texture) only normal mass hierarchy is possible (texture A). Results are given in Fig. 
In this figure, the second row gives solutions when experimental errors are shrinking two times, without changing 
central values a:^*^"'. We can see, that regions decreased drastically. It means that many of the solutions found are 
placed nearby their central values and many of peripheral solutions connected with edged errors are cut away by our 
rejection rule given by the a parameter. In Fig.|B|the solutions are transformed to the physical parameters connected 
with mixing angles. 

Coming from the numerical solutions of the mixing matrix U into the language of mixing angles, the standard form 
of the neutrino mixing matrix represented by three angles 6*13, ^23 and 9i2 has been used (see e.g. [5|, CP phase 
(5 = 0). If a specific hierarchy of the neutrino masses is chosen, the mixing angles defined by oscillation probabilities 
can be restricted to the first quadrant [22 ■ However, in Fig. (and also Figs. I8I1U|I we consider general approach with 
unbounded mixing angles. We can see that some of histograms are not symmetric against sin^^ — 0. In Fig. 6 we 
keep the plot of sin023, though for the experimentally and phenomenologically interesting case a — Q the neutrinoless 
double beta decay is independent of ^23 (see e.g. 24, 25] ). It means that the same histogram results for the case 
when a 7^ 0. Values of suiOij parameters in Figs. I6I8I1QI correspond to the 3 a confidence level. We can see in Fig. El 
that 013 around zero is preferred for this texture. Histograms for all other cases with a one zero parameter can be 
found in addition in ^6]. There are analytical considerations for the case a = in the literature, see e.g. |2J,|23. It 
is easy to prove that in the case a = 0, and for Oij, = a simple relation holds 



mi 



m2 



= tan^6'i2. (5) 



We have checked that our numerical solutions fulfill this relation. 

Let us discuss results for the case where two zero parameters are presented in Eq. |21 

If n of six independent parameters in Eq.|31 vanish then one can find „i(g^l„S| possible textures. For n — 2, we have 
15 cases. However, we have found that, to be in agreement with present experimental data (i.e. a — \, ia c.L), seven 
cases are possible. All possible two zero textures are shown in Table E] In the case of the normal mass hierarchy, 
we have got the so called Ai and A2 textures |23| with a, 6 = and a^c — 0, respectively. B textures have both 
normal and inverse hierarchies and C texture has inverted one (see Table llTT|) . From the general case (inverted mass 
hierarchy in the last section) we have concluded that d, e, / parameters can be zero, however, it appears that only 
situation when d = and / = is possible, all other combinations, e.g. d = and e = are excluded (to come 
to this conclusion, altogether around 10^ randomly generated points have been used for each degree of freedom). In 
Table HhI we recalculated solutions for M^ to get physical values of neutrino masses. The last column shows results for 
the ao parameter for which no solutions have been found. As already discussed, smaller ag, more parameters around 
zero are allowed. The calculation is stopped when 10^ points for each degree of freedom for the Monte Carlo random 
generator is not able to find solutions in agreement with Eq. 0| The one before last column in the Table IIIII shows a 
prediction for the neutrinoless double beta decay parameter < mp poi, >. We will consider this important parameter 
in more details in a separate paper, where also possible CP violating phases will be examined. 
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TABLE IL The table shows all two zero textures which are in agreement with current experimental data. 

From the Table Ulll we can deduce (see the best values in the MASS MEAN column) that neutrino masses fulfill 
relations: 
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TABLE III: Masses and effective neutrinoless double beta decay mass parameter < mfjfjQu > for allowed textures with two 
zeroes and one with one zero. The last column shows the parameter qq for which schemes have no positive solutions. 
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TABLE IV: This table shows a;^^"' values obtained from numerical solutions for two zero textures. It appears, that cases Ai 
and A2 coincide with solutions for one zero texture with a — 0. 



mi < 7712 << 7713, texture A (6) 

7TT-3 << mi < 7712, texture C, (7) 

in agreement with analytical considerations given in [l^l and [lj|- However, for textures B the situation is more 
complicated. In the above cited papers a degeneration of neutrino masses for textures B have been inferred. However, 
as we can see from Tabic iHll for textures B not only the complete degeneration of neutrino masses 77ii ~ 7712 ~ m^ is 
allowed. Actually, both inverted and normal mass hierarchies are possible. We think that the reason for discrepancies 
lies at assumptions or approximations made in the analytical approach which cut off some still possible solutions. 
Let us note, that most of analytical results to which we refer include cases with CP phases. Here we have restricted 
ourselves to the real neutrino mass matrix cases, so we can expect even wider spectrum of solutions when complex 
mass matrices will be included. 

In Table HVI the results with the smallest x^ values found by our adaptive Monte Carlo procedure are written. We 
can see that some of them "walked away" from corresponding values in Table ^ 

Finally, in Fig.[7|we visualize results for the C texture. Corresponding histograms are shown in Fig. |S1 In Fig. |3 
and Fig. ^] the histograms for A and B textures are given, respectively. We can see that for schemes A, nonzero ^13 
angles are favorable. Similarly to the case with a one zero texture, additional set of figures can be found in |2y|. 

We have found, that there are no allowed textures when 3 (and more) parameters of the matrix Mi, are exactly 
zero. This situation can be different when CP violating phases are taken into account (work in progress). 

IV. CONCLUSIONS 

We have analyzed the special case of Majorana neutrino textures which give CP conserving effects. Regions of 
parameters have been found which are in agreement with present experimental data. Allowing for random generation 
of parameters, histograms suggest some preferable regions of parameters. Moreover, the results show that neutrino 
textures with some zero elements are the most probably solutions. 

Let us summarize the most important points: 



• 



• 



for the general case, some elements of the neutrino mass matrix M^ are likely to be around zero, 
there are no possible numerical solutions for M^, textures with number of zeros 71 > 3 (at 3cr c.L), 



• there are seven two zero textures which give results in agreement with present experimental data, some of them 
can have both normal/degenerate and inverse/degenerate mass hierarchies, some of them have only normal, or 
only inverted mass hierarchies, 

• textures B give small values of sin0i3, 

• cases with mppQn = have got only normal mass hierarchy and they all imply similar results. 
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APPENDIX A: DESCRIPTION OF THE ADAPTIVE MONTE CARLO ALGORITHM 

The main purpose of the algorithm is to find values of parameters of neutrino mass matrix Af^, Eq. |21 which 
lead to results which are within ranges of current neutrino experimental data. It is obvious that it is not pos- 
sible to check numerically the all possible combinations of input parameters, the Adaptive Monte Carlo (AMC) 
is the best solution. The program works in two steps. Firstly, it generates randomly input parameters and 
verifies the obtained output with experimental data. Thus the basic set of parameters is obtained. Such sets 
are represented as points on plots, see e.g. the left plot in Fig. ^ The user can specify how many points he 
wants to obtain. For a typical run it is Ngcat = 100- As soon as the first step ends, AMC algorithm begins. 
It uses points which are obtained during the first step and searches for a solution with the lowest value of y^ 
around every of these points. Second step consists of Nu iterations (enlargements), and for every iteration a 



diagonalization of the neutrino mass matrix is performed by generating randomly points Ncaic times. As a result 
of the second step, additional points are obtained which form regions of allowed parameter space, e.g. the right 
plot in Fig. nj One can notice that these regions look denser (more points) in the center than near the edges. It 
is a result of an algorithm which zooms into the points with the lowest values of x^- Below is a sketch of the algorithm: 

1. First step: Scattering 

• Random generation of input parameters 

For the nth iteration (71 = 1,2,...) the program generates set of input parameters: neutrino mass matrix 
elements. Algorithm stops the loop when n = Nscat- Because in this work we consider Majorana neutrinos 
then neutrino mass matrix is symmetric so this implies that diagonal elements of the mass matrix are real. 

• Diagonalization of the neutrino mass matrix M^, 

With mass matrix input parameters chosen, the algorithm finds eigenvectors (neutrino masses) and nor- 
malized eigenstates (columns of the neutrino mixing matrix) 

• Comparison with experimental results and saving allowed parameters 

The program compares the obtained values with experimental results (neutrino mass squared differences 
and absolute values of neutrino mixing matrix elements) with a help of the Xi function 

(cent _ „ .\2 

where Xi is a calculated value and a;^*^"*, a]^ are experimental best fit and uncertainty values, respectively. 
The constant a is a parameter which controls decreasing of an assumed error, Eq. 0] If calculated x^ is 
smaller than the demanding number, the parameters (points) are saved. 

The purpose of this step is to scatter randomly points defined by the neutrino mass matrix parameters in a 
hyperspace volume where a number of dimensions is equal to the number of input parameters. 

2. Second step: The Adaptive Monte Carlo 

• Reading obtained points 

For the n-th iteration (n = 1, 2, . . . , A'scat) program reads set of parameters (points) obtained in the 
previous step and sets their values as the central values x^*^"'. For the i-th iteration [it = 1,2,..., iV^t), 
the algorithm sets range of input parameters (which will be randomly generated in the next step) : 

X^''±^^t5^. (A2) 

where i numerates input parameters, 5^ is an initial range and (^a is a function defined as follows: 

The constant number 0.6 and the function \/it are set empirically. Similarly as in the first step, the 
program executes the following (a)-(c) tasks N^aic number of times: 

(a) Random generation of input parameters 

(b) Diagonalization 

(c) Comparison with experimental data and saving successive cases 

• Setting new central values 

When tasks a) — c) are finished, the algorithm chooses these set of input parameters for which the calculated 
value of x^ is the lowest: 

X^ = Ex?- (A4) 

z=l 

and sets it as "the best set" , i.e. as a new central value x^^"*. Next, the program repeats above calculations 
with the new value of it which is it + 1. Note that while increasing it, the value of function ^it decreases, 
so the range in which we randomly generate input parameters will be smaller and smaller as it increases. 
Also, since every time the best x^ is chosen, the solution moves towards the most probable solution for a 
considered set of experimental results. 
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FIG. 1: Three dimensional plots of allowed parameters found by the AMC procedure. On the left plot there are points obtained 
firstly by generating random parameters which fullfill relation Eq. lAll on the right plot the points are denser as AMC looks 
for additional solutions in a vicinity of the parameters obtained in the first step Eas lA2IA4l These plots show solutions for the 
Al texture. 
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FIG. 2: Allowed regions of parameters for the neutrino mass matrix M^ with present experimental data (TableQ] Sa level), the 
general case with normal mass hierarchy (upper plots) and inverted mass hierarchy (bottom plots). 
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FIG. 3: Frequency spectrum for the elements of the neutrino mass matrix M^: the general case with normal mass hierarchy. 
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FIG. 4: Frequency spectrum for the elements of the neutrino mass matrix M^: the general case with inverted mass hierarchy. 
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FIG. 5: Allowed regions for the mass matrix with a = (A texture). The first row shows plots with a — 1 (present data, So- 
level), the second row shows results for a = 2. 
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FIG. 6: Histograms of neutrino rotation angles sinSi2, sinS23 and sin 6*13 for neutrino mass texture with one zero a — 0. The 
histogram for sin ^23 does not depend on a and is the same for a 7^ 0. 
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FIG. 7: Allowed regions for the mass matrix with d, / = (C texture). The first row shows plots with a — 1 (present data, So- 
level). The second row shows results for a = 2. 



leooo 




i — 1 — ' 


1 ' , ' 


— 1 — ' 


1 1 1 


' — n 




1 — ' — ' 


1 ' 


1 i r- 


— 1 — 


' 1 






6000 




-1 ^-^ • 1 ' 1 n 


14000 












1 


20000 




















_^^M 


12000 












1 


16000 
















4000 




JUL 


10000 
6000 












I 


12000 
















3000 


i 
















L 


6000 












k 


8CO0 
















2000 


_l 
















■ ' 


4000 

2000 














1 


4nnn 




L 








, 






1000 



jff 




















k 


o.e 


-0,4 


-0,2 0.0 


0,2 


Q,4 


0,6 


-0.6 


-0,3 


0,0 


0,3 


0,6 





9 


-0.14 


),07 0,30 07 


0.14 








Bine,, 














sine„ 


















S 


.in 


6 













FIG. 8: Histograms of neutrino rotation angles sin 6*12, sin&23 and sin 6^13 for neutrino mass texture with d, / = (texture C). 
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FIG. 9: Histograms of neutrino rotation angles sin 612, sinS23 and sin 6^13 for neutrino mass texture with a,b = and o,c = 0. 
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FIG. 10: Histograms of neutrino rotation angles sinSi2, sin 623 and sin^is for neutrino mass textures B (all cases) 



